Spectroscopy of SMC Wolf-Rayet Stars Suggests that 
Wind- Clumping does not Depend on Ambient Metallicity 
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ABSTRACT 

The mass-loss rates of hot, massive, luminous stars are considered a decisive 
parameter in shaping the evolutionary tracks of such stars and influencing the 
interstellar medium on galactic scales. The small-scale structures (clumps) om- 
nipresent in such winds may reduce empirical estimates of mass-loss rates by an 
evolutionarily significant factor of > 3. So far, there has been no direct observa- 
tional evidence that wind-clumping may persist at the same level in environments 
with a low ambient metallicity, where the wind-driving opacity is reduced. Here 
we report the results of time-resolved spectroscopy of three presumably single 
Population I Wolf-Rayet stars in the Small Magellanic Cloud, where the ambient 
metallicity is ~ 1/5Z & . We detect numerous small-scale emission peaks moving 
outwards in the accelerating parts of the stellar winds. The general properties 
of the moving features, such as their velocity dispersions, emissivities and aver- 
age accelerations, closely match the corresponding characteristics of small-scale 
inhomogeneities in the winds of Galactic Wolf-Rayet stars. 
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1. Introduction 

Wolf-Rayet (W-R) stars, as a final evolutionary phase of hot, massive stars [Mmit > 
25M for aZ~ Z & ambient metallicity), are especially sensitive to all phenomena which may 
influence their mass-loss rates. The fast, dense winds of W-R stars are driven by radiation 
pressure exerted on multiple lines of mainly heavy elements. Hence, W-R mass-loss rates 
should be sensitive to the ambient metallicity content, as well as to the locally enhanced 
chemical composition of the wind. The former sensitivity, though suggested from general 
principles (Lamers & Cassinelli 1999), has escaped detection until recently (see Crowther 
2006, Grafener & Hamann 2006, and references therein). 

Line-driven winds are inherently unstable (Lucy & Solomon 1970), being constantly 
fragmented by numerous embedded shocks (Dessart & Owocki 2005). Such fragmentation 
may change the local ionization balance and create a local non-monotonic velocity field, 
thus providing a strong feedback to the driving force. In the past, wind diagnostics relied 
on models of smooth, homogeneous winds, until numerous spectroscopic observations of 
Galactic W-R and OB stars (Moffat et al. 1988; Robert 1992; Lepine & Moffat 1999; 
Eversberg et al. 1998; Bouret et al. 2005; Fullerton et al. 2006 and references therein) 
demonstrated the omnipresence of wind-embedded clumps, usually taking on the form of 
discrete density enhancements outmoving with an accelerating wind. As an immediate and 
long-lasting impact, the relatively simplistic treatment of structures in the otherwise highly 
elaborate models of hot-star winds (Hillier & Miller, 1999; Puis et al. 2005; Hamann et 
al. 2006 ) resulted in a consensual downward revision of mass-loss rates by a factor of > 3, 
confirmed by numerous independent observations (e.g., Moffat & Robert 1994). 

Though wind-clumping was widely anticipated to operate also in low-Z environments 
(Hamann & Koesterke 1998; Crowther et al. 2002; Bouret et al. 2003), any direct proof was 
lacking. Here we report the detection of outmoving spectral features in the line profiles of 
3 SMC W-R stars. We compare the general characteristics of these clumps to those in a 
Galactic sample of W-R stars, finding them to be strikingly similar. 

2. Observations 

We targeted three presumably single (Foellmi et al., 2003) early-type WN stars in the 
SMC: SMC-WR 1 (WN3ha), WR 2 (WN5ha) and WR 4 (WN6h), thus forming a represen- 
tative sample of the small SMC W-R population (12 known Population I W-R stars, most 
of them early-type WNs: Massey et al. 2003). We monitored the stars in continuous, ~lh 
40 min-long loops for two consecutive nights in August, 2006 (between HJD 2,453,974.538 
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and ...75.877), alternating among the 3 targets: WR 4 - WR 2 - WR 1... . We used 
the UVES spectrograph at the ESO-VLT-UT2 8-m telescope (Kueyen), sampling the re- 
gion A3927 — 603lA . During the routine reduction we experimented with different binning 
factors, finding a reasonable compromise between desired spectral resolution, < 0.5A /pix, 
and signal-to-noise ratio, S/N > 100, by binning the available spectra to 0.25A (the region 
around Hell 4686A ) and 0.3lA (Hell 5412A ) bins. With 40, 30 and 20 min exposures, this 
resulted in 5-6 spectra/ night /star, with comparable signal-to-noise ratios S/N ~ 120, ~ 140 
and S ~ 160 (±20) for WR1, WR2 and WR4, respectively, estimated from the adjacent 
continua around measured spectral lines. 



3. Properties of the Clumps 

The high quality of the spectra allowed us to apply a rather straightforward procedure 
to detect the line-profile variability. We rectified the profiles of all relatively prominent Hell 
lines (the only lines strong enough to provide sufficient precision, namely, Hell 4686, 4859 
and 5412 A ) by linear interpolation between the line-free regions located in the immediate 
vicinity of a given line. Then we created unweighted average profiles of each line for each 
night of observation and subtracted these averages from the individual profiles. We found 
that all the major emission lines in all 3 stars can be considered as variable at a level 
substantially exceeding the predicted instrumental noise level. 

For a proper comparison of the new results with previous studies of inhomogeneities 
in the winds of Galactic W-R stars, we closely follow the already developed approaches 
specifically designed to reveal basic properties of clumps. 

(i) First of all, we attempt to detect similarities in the variability patterns of major 
emission lines by cross-correlating the difference (individual - average) spectra of different 
line transitions. We use the average cross-correlation functions (see Lepine et al. 2000, 
equations 7, 8): namely, we produce an unweighted average of the individual (transition-to- 
transition) cross-correlation functions for a given star and a given spectral window (usually, 
±t>oo, centered on the transition's Ao + AA, where AA corresponds to the SMC's systemic 
velocity). The maxima of the average cross-correlation functions reach 0.2-0.4, which is 
comparable to I max = 0.2 — 0.8 measured in the CIII and CIV profiles of the Galactic 
star WR 135. Taking into account that (a) the considered transitions arise from different, 
though overlapping, line-forming regions, (b) there is some hydrogen present in the A4860A 
blend, and (c) there might be some transition-dependent optical-depth effects (see below), 
we did not anticipate a very high degree of correlation. Nevertheless, for general consistency, 
plus a slight gain in S/N, we combine the appropriately shifted (to the rest-frame of the 
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Hell 4686 transition) difference spectra of the three lines, smooth them with a gaussian 
filter (FWHM=7 pix), linearly time-interpolate the neighboring spectra in order to fill the 
unavoidable gaps, and plot the results in Fig. 1. The continuous, V-shaped structures present 
in the difference spectra of the SMC W-R stars are remarkably similar to the variations seen 
in all appropriately studied Galactic WR stars (Lepine & Moffat 1999), with the latter 
firmly linked to the presence of numerous small-scale inhomogeneities moving outwards in 
the accelerating parts of the stellar winds. 

(ii) We proceed with measurements of the clump properties closely following the ap- 
proaches developed in Lepine & Moffat (1999) and in Lepine et al. (2000). Namely, we 
use the degradation function to estimate the average acceleration of the outmoving emission 
features (Fig. 2) and the net intrinsic variability levels, a intr , in different transitions (Fig. 

3) . We remind the reader that the degradation function (Lepine et al. 2000, equation 2) 
measures the mean standard deviation between spectra separated by a [variable] interval At. 
In addition, this function also accounts for a possible motion of the clumps (predominantly, 
away from the line's center), by including the average acceleration as a free parameter. The 
levels of net intrinsic variability are estimated across emission-line profiles (Lepine et al. 
2000, equations 4,5) by using the variability in line-free continuum regions adjacent to emis- 
sion lines and treating it as a source of instrumental noise, then adjusting this noise level by 
allowing for the difference in count statistics across the profile of a prominent emission line. 

(iii) In order to compare our measurements of the individual fluxes and full-widths at 
half-maxima of the detected clumps with corresponding values for Galactic WR stars, we 
follow the general approach and published measurements from Robert (1992), thus using only 
the Hell 5412A line. Namely, we measure the fluxes, maximum intensities and FWHMs from 
the least-square fits of Gaussian profiles to the recognizable features in the difference-spectra 
(individual spectrum - night average). Fig. 2 shows that, within the statistical uncertainties, 
which are larger for the SMC stars mainly due to the relatively lower S/N, the properties of 
the Galactic clumps are indistinguishable from those in the SMC sample. Note that for the 
estimates of the clump properties we use the unsmoothed difference-spectra. 

The broad spectral coverage of the data allows us to address the question of optical 
depth of the clumps by comparing the intrinsic variability patterns in different transitions. 
One should note (Fig. 3) the most instructive contrast between Hell 4686A (transition 3- 

4) and Hell 5412A (transition 4-7). In order to make quantitative inter-comparisons of the 
intrinsic variability levels, we normalise <Ji ntr by the line-profile intensity of the corresponding 
transition. Both SMC-WR 1 and WR 2 show a reasonably consistent behavior in both 
transitions (Fig. 3; compare to Fig. 6 from Lepine et al. 2000), if one takes into account 
the pronounced difference in the extensions of the line-forming regions: compare the velocity 
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extensions of the average profiles of Hell 4686A to their Hell 4859 and 5412A counterparts. 
Surprisingly, the variability patterns of WR 4 come as a sharp contrast to the predictably 
congruent behavior of the patterns in WR 1 and WR 2. 

4. Discussion and Conclusions 

First of all, we concentrate on the intrinsic line variability as related to optical-depth 
effects introduced by the clumps. Relying on the differences in the total optical depths of the 
studied transitions and assuming that these are also applicable to the population of discrete 
density enhancements, one expects to detect different cr intr dependencies across the profiles. 
Most notably, the significant optical depth of Hell 4686A may result in a blue-shifted peak 
of a intr relative to Hell 4859A and Hell 5412A . Quite predictably, SMC-WR 1 and WR 2 
show little, if any, line-opacity effects: note the pronounced (Fig. 3) absorption components 
in the profiles of Hell 4859A and Hell 5412A , which are generally linked to low optical 
depths of the winds. Concurrently, Lepine & Moffat (1999) find only mild clump-related 
optical-depth effects in the [optically thin] CIII 5696A line formed in the winds of Galactic 
WC stars, along with the possible presence of opacity-related effects in Hell 5412A seen in 
the Galactic WN stars. Hence, for WR 1 and WR 2 one may conclude that either (a) the 
clumps are optically thick in all transitions and the winds have rather low volume filling 
factors (e.g., Owocki & Cohen 2006), or (b) all clumps are optically thin. Unfortunately, our 
present limited data set does not allow one to choose between these two alternatives. The 
only noticeable difference between WR 1 and WR 2 is the secondary maximum in ai ntr of 
Hell 4859A for WR 2, which may be caused by the presence of relatively high quantities of 
hydrogen in the wind. 

WR 4 tells a different story. The intensity of the Hell 4686A transition, as well as the 
lack of developed P Cygni absorptions, calls for a denser wind. If the global properties of 
the wind are related to the properties of clumps, then one may detect different wavelength 
dependencies of a intr . This clearly pertains to WR 4, where the blueward displacement of 
the Hell 4686A peak may point to a high optical density of clumps. However, the disparity 
of the Hell 4859 and 5412A distributions prompts an alternative explanation: the spectral 
variability could be, in part, related to the presence of co-rotating interacting regions (CIRs: 
Massa et al. 1995; Morel et al. 1999, and references therein), as seen in several Galactic W-R 
stars. One may notice the synchronized redward migration of weak emission features in the 
first half of Night 1 (Fig. 1, dashed lines: A ~ 4673 shifting to A ~ 4676, ~ 4689 gradually 
moving to ~ 4693, and ~ 4698 moving to ~ 4702), as well as the gradual disappearance of 
the broad emission feature at A ~ 4690 during Night 2, both phenomena reminiscent of a 
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CIR. In addition, along with general expectations, WR 4 is a periodic photometric variable 
(Foellmi et al. 2003: P=6.55 d), and not an obvious binary. This period is then likely related 
to rotation. However, this does not mean that WR4 lacks clumps. 

One should also notice the lack of any dependence of clump FWHMs on spectral class 
or terminal velocity of the W-R winds (Fig. 2, lower section). The recent comprehensive 
analysis of turbulence in colliding supersonic flows (Folini & Walder 2006) shows that the 
velocity dispersions in the shock-bound cold dense layers (the clumps?) depend only on 
the upwind Mach-numbers, M u . We consider a 'generic' W-R wind with = 2000 km/s, 
estimating the dispersion of turbulent velocities as cr(vt U rb) ~ (0.1— 0.2)foo (Prinja et al. 1990; 
Marchenko & Moffat 1999). We assume an 'asymptotic' wind temperature T = 2- 10 4 K. This 
provides an estimate for the sound speed in the clump-formation zone, ~ 5 — 100-R* (see the 
models in Lepine et al. 2000). Using M rms ~ rf l l 2 M u with density contrast 77 ~ 30 (Folini & 
Walder 2006), one may estimate the typical width of a clump as: FWHM ~ (1.4 — 2.7) A , or 
90-170 km/s at A = 4686A . This rough estimate falls tantalizingly close to the measurements 
from Fig. 3, which give FWHM = 100 — 300 km/s. One may match these figures even 
more closely under the assumption that the hydrodynamically-modeled density contrasts fall 
below the pre-specified value of rj = 30, by adopting rj ~ 10 as frequently used in models of 
W-R spectra (Hamann & Koesterke 1998). 

Hence, quite justifiably, we favor the notion in which the clumps have the ~same internal 
velocity dispersion defined by compressible turbulence in a supersonic flow, but different 
emissivity volumes (cf. Lepine & Moffat 1999). 

Concluding our study, we can state that: 

(a) despite the differences in ambient metallicity, both the winds of Galactic and SMC W-R 
stars show a clear presence of small-scale structures; 

(b) the general properties of these structures, namely average acceleration rate, average indi- 
vidual flux and average velocity dispersion inside the clump, are similar for both populations 
of W-R stars; 

(c) point (b) comes as a big surprise, considering that the wind-driving force (and thus the 
radiatively-driven instability) sensitively depends on the presence of heavy elements (Vink 
& de Koter 2005); 

(d) a wind-clumping factor must be included in the models of massive-star winds in low-Z 
environments, just as in high-Z environments. Our understanding of the first population 
of stars in the early Universe, mostly all very massive, as well as processes of enrichment 
of the interstellar medium, will depend on reliable estimates of mass-loss rates, i.e., the as- 
sumed wind-clumping factors. Indeed, wind-clumping could be rather extreme in a low-Z 
environment: one may recall f=0.01 used in the models of MPG 324 (04 V((f)) star in 
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SMC: Bouret et al. 2003), and f=0.06 for Br43 (WC4 star in LMC: Hamann & Koesterke 
1998), to be contrasted with the more modest / ~ 0.1 — 0.3 values for Galactic WR stars. 
However, note that the recent data on Galactic OB stars also provide / <0.1 (Bouret et al. 
2005; Fullerton et al. 2006). This may signify that: either the filling factors depend on the 
evolutionary status of massive stars, or the invoked wind-clumping factors are biased by the 
model approach and the choice of diagnostic lines. The latter, in turn, may be related to a 
[strong] radial dependence of the wind-clumping factor. This obviously calls for a thorough 
investigation which is far beyond the scope of this paper. 

The extreme clumping may dramatically decrease estimates of the mass loss rates, thus 
helping to retain the expected high rotational velocities for the progenitors of the long/soft 
7-ray bursts, as well as affect the appearance of their immediate environments, thus the 
observable reaction (echo) on the incoming burst. 

This paper is based on observations made with ESO Telescopes at the Paranal Ob- 
servatory under programme ID 077.D-0029. C.F. thanks C. Esparza for the efficient and 
friendly technical support during observations. AFJM is grateful for financial assistance to 
NSERC (Canada) and FQRNT (Quebec). FM acknowledges support from the Alexander 
von Humboldt foundation. 
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Fig. 1. — Grayscale plots of time- interpolated and smoothed difference (individual - night- 
average) spectra of WR1 (left panels), WR2 (central panels) and WR4 (right panels) for 
Night 1 (lower panels) and Night 2 (upper panels). The intensity ranges are -0.03 (black) to 
+0.03 (white) in the local continuum units (= 1). The dashed lines in the lower-right panel 
mark the features presumably related to CIRs (see text). 
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Fig. 2. — The upper panel shows the average accelerations of outmoving clumps in the winds 
for SMC stars (filled circles), as well as the measured ranges for Galactic WC (open squares) 
and WN (open triangles) stars (see Lepine & Moffat 1999). The data are arranged by the 
corresponding terminal velocities of the WR winds. The middle panel shows the average 
fluxes (F) of the clumps detected in the Hell 5412A line and normalized by the maximum 
intensity of this emission profile (I max — 1), in order to be compared with the measurements 
of Robert (1992). The lower panel shows the average FWHMs of the clumps observed in the 
Hell 5412A profile. Note that WR 134 was omitted from the lower panel due to the very 
high FWHM ~ 800km/ s (see explanations in Lepine & Moffat 1999). The estimates 
are rather uncertain for WR 1 and WR 2 (P. Crowther, priv. comm.); for WR 4 they come 
from Crowther (2000) and Willis et al. (2004). 
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Fig. 3. — The intrinsic variability across the line profiles in WR 1, WR2 and WR 4 (lower 
panels, left-to-right): full lines correspond to the Hell 4686A line, dotted lines show Hell+H 
4860A , dashed lines trace Hell 5412A . The upper panels show average rectified profiles of 
the transitions. 



